Extended Moment Formation and Second Neighbor Coupling in Li2Cu02 
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Comprised of ferromagnetic edge-sharing Cu02 chains that order in antialigned fashion at Tjv=9 
K, Li2Cu02 is found from local spin density calculations to display several surprising characteristics: 
(1) the ordered moment/f.u. of 0.92 is the largest for any low dimensional cuprate system, in 
agreement with experiment, (2) 40% of this moment lies on the neighboring O ions, making it the 
largest oxygen moment yet reported, and (3) the second neighbors couplings are larger than nearest 
neighbors couplings. All of these phenomena arise naturally due to a well defined eff'ective dy^ 
type orbital that includes very strong O pa character. We interpret the large moment as surviving 
reduction by quantum fiuctuations due to extension caused by the d ~ p hybridization. 



Low dimensional (ID or 2D) spin S=i materials are 
critical systems displaying a variety of phenomena in 
crystals. These range from spin gap behavior in CaV409 
101 to anomalous spin physics in the one-dimensional spin 
chains and ladders. |^ Many of these characteristics de- 
pend crucially on specific structural or chemical bonding 
features. |^ The magnetic coupling is particularly sensi- 
tive, with simple nearest neighbor (nn) exchange vary- 
ing from large and antiferromagnetic (AF) to small and 
ferromagnetic (FM) when the metal-oxygen-metal angle 
(j) varies from 180° to 90°. Weak intrachain couplings 
and competing exchange couplings can lead to frustra- 
tion, magnetic ordering, spin gap behavior or spin-Peierls 
phase formation. 

Copper oxides play a very important role due to the 
various possibilities of linking their fundamental unit, a 
(often slightly distorted) CUO4 square. Three general 
types of arrangement can be found in these systems, clas- 
sified in terms of the oxygen squares sharing corners (as in 
the high Tc planar compounds and Sr2Cu03), edges (as 
in CuGeOa, La6Ca8Cu2404i, and Li2Cu02), 0, or both 
(as in SrCu02). In edge sharing CuGeOa, for example, 
there is moderately strong antiferromagnetic (AF) nn 
couphng (J«150 K), [| but a spin-Peierls transition oc- 
curs only at 14 K. The typical example of corner-sharing 
chains is Sr2Cu03, where the system orders antiferro- 
magnetically Q with a low transition temperature of 5 
K and a small induced magnetic moment (0.06 /is) in 
spite of very large interactions between ions (J«2200 K). 
Despite considerable progress a clear understanding of 
the magnetic behavior of the Cu'^^ ion in several regimes 
is still lacking. 

Here we report a new aspect of spin behavior in edge- 
sharing systems revealed by spin-polarized local den- 
sity approximation (LDA) studies of a ID system 
Li2Cu02. In this compound neutron scattering Q in- 
dicates three dimensional AF ordering at 9 K arising 
from the antialignment of FM chains. The experimen- 
tal moment of 0.9 /is per cell was attributed completely 
to the Cu ions. Based upon experience in the undoped 



two-dimensional cuprates {viz. La2Cu04) where LDA is 
unable to obtain any moment whatsoever on the Cu ion 
§ , it might seem that LDA is unlikely to produce a mag- 
netic Cu ion or an insulating system. Due to this problem 
with corner-sharing Cu-0 planes, very few spin-polarized 
calculations on cuprate compounds have been reported. 
However, we find that LDA predicts Cu in Li2Cu02 to 
be robustly magnetic, allowing us to obtain the relative 
energies and electronic properties of the system with FM 
chain with both AF and FM coupling between chains, 
the AF chain, and the unpolarized (PM) system taken 
as reference. The oxygen ions play a fundamental role in 
the band dispersion and the magnetism, and carry a mag- 
netic moment approaching 0.2 ^lb per atom, the largest 
O moment yet reported. Previous reports of O moments 
lie in the 0.02-0.10 ^ib range. ^ 

Li2Cu02 is orthorhombic and belongs to the cate- 
gory of edge-sharing compounds, with one-dimensional 
Cu02 ribbons carrying the Cu chains along the h axis, 
arrayed in a body-centered fashion in the a — c plane 
(Fig. p. The Cu-O-Cu angle = 94°, intermediate be- 
tween the two other edge-share compounds with very dif- 
ferent characteristics: GeCuOa (spin-Peierls, = 99°) 
and La6Ca8Sr2404i (FM, (f) — 91°). Distances between 
two Cu ions, 2.86 A along the chain, 3.65 A in the a — b 
plane and 5.23 A in the diagonal direction, do not reflect 
the relative coupling strengths, as we explain below. 

Calculations were done using the linearized augmented 
plane wave (LAPW) method, jl^ which makes no shape 
approximations for the density or potential. The sphere 
radii used in fixing the LAPW basis were chosen to be 
2.00 a.u for Cu and Li and 1.65 for O. Local orbitals (Cu 
3p; Li Is; O 2s) were added to the basis set for extra flex- 
ibility and to allow semicore states to be treated within 
the same energy window as the band states. The plane 
wave cutoff corresponded to energy of 23.5 Ry resulting 
in 640 LAPWs per formula unit. Self-consistency was 
carried out on k-points meshes of 512 points in the Bril- 
louin zone for the compounds which need a single unit 
cell calculation and 256 points when we considered an 
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FIG. 1. Crystal structure of Li2Cu02, for the AF unit cell 
containing two Cu02 ribbons. 

AF arrangement for the chains. 

The paramagnetic system has an odd number of elec- 
trons per unit cell and thus is metallic. There is a single 
band in the range of 1 eV around the Fermi level split 
off from the rest of the p — d band complex, similar to 
what was found in the other Cu02 edge-sharing com- 
pounds CuGeOa fill and NaCu02. This isolated 
band shows up as two bands in Fig. b|, where a doubled 
cell with two chains has been used for comparison with 
the AF bands (see below). The analysis of the partial 
density of states shows that only Cu dj^2-0 Pa (the Py±pz 
combination directed towards the Cu site) are present in 
this band. 

The geometry of this Cu02 edge-sharing chain leads 
to a simple description of the important band, which we 
expect (and find) to involve an antibonding combination 
of Cu d and O p orbitals. The atomic orbital basis in 
a primitive cell can be chosen as the Cu d orbitals, the 
CT-type O Pfj orbitals on each of the Jour neighboring O 
atoms which strongly overlap the dy^ orbital, and half of 
the out-of-plane Px orbitals that are non-bonding and lie 
well below Yip (as do all d orbitals except dyz ) . The Pir or- 
bitals in the y — z plane of the ribbon are Pa with respect 
to a neighboring Cu and belong to the next unit cell. 
The dyz and the four pa orbitals can be decomposed into 
five hybridized combinations: one bonding combination 
"Dyz and one antibonding combination I?*^ of dyz and the 
combination of the four p„ orbitals with dyz symmetry, 
and three other po.-only combinations of lower symme- 
try. Vyz and V*^ are split strongly, leaving V*^ at the 
Fermi level and T^yz 5 eV below, as shown in the den- 
sity of states plot in Fig. ||. The others lie around -4 eV 
below Yip and are not of interest. The general behav- 
ior of the coupled dyz — Pa- cluster can be modelled with 
Ed — — 1.5,ep = —4, {dpa) = ±1.15, [ppa) = 0.25, (all in 
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FIG. 2. Band structure of a) paramagnetic Li2Cu02, b) 
a zoom of the region around the Fermi energy, and c) 
for antiferromagnetic alignment of FM chains. F, X, Y 
and Z points correspond to the standard denomination, 
A = (tt, tt, tt), R = (tt, 0, tt) and T = (0, tt, tt); all in the dou- 
ble cell Brillouin zone. 

eV) . For these parameters the V^y density is 70% on the 
Cu and 30% on the four O ions. 

The active orbital V*^ shown (schematically) as \'Dy^\'^ 
on next nearest neighbors in Fig.^, is an effective d^^-type 
orbital centered on each Cu ion but extending strongly 
to the neighboring O sites. Symmetry allows direct "DDtt 
overlap, and therefore hopping amplitude along the 
chain. Due to its parentage, however, it is clear that 
the main contribution to the overlap arises from the O 
ion region. If the O quadrilateral were perfectly square 
{(f) — 90°) the O Pa orbitals directed toward the two 
neighboring Cu ions would be precisely Py i Pz- These 
combinations are orthogonal, so reduces to direct d- 
d overlap and will be very small. When the Cu-O-Cu 
angle is not exactly 90°, the Pa orbitals are no longer 
orthogonal and the overlap (and tj^) increases. 

The effective Hamiltonian therefore reduces to a single 
orbital (T)*^) per cell. The dispersion of the V*^ band 
cannot be fit simply by nn hopping = ti along the 
chain and t[ between neighboring chains (in roughly the 
X ± z directions), but requires as well both next nearest 
neighbor (nnn) hopping terms t2 and t'2- Consideration 
of the orbitals on second neighbors indicates why 
this is so: Cu-O-O-Cu coupling along the ribbon becomes 
important because of 0-0 coupling and because the nn 
hopping is so small, and similarly for interchain hopping 
along the diagonal. The manner of second neighbor over- 
lap along the chain is clear in Fig.^. The values 

nn : <i = —63 meV , t'l = — 16 meV, 
nnn : <2 = ^94 meV , ^2 ^ 44 meV, 
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FIG. 3. Density of states for the paramagnetic system and 
for AF alignment of the chains. In the PM case the partial 
density of states for the Cu dyz is shown as the shadowed 
region, illustrating the 5 eV splitting of the "Dyz and 
bands. 

(note that second neighbors values are larger than near- 
est neighbors values), provide an excellent fit to the dis- 
persion except near the (7r,7r, k^) line (not shown) where 
some additional interaction would be required to obtain 
a perfect fit. 

This emergence of an effective single-band system with 
a simple, strongly hybridized orbital is one of the cleanest 
in Cu-0 systems. It implies that the (ferro- or antiferro- 
) magnetism of this compound should be interpreted in 
terms of this extended V*^ orbital instead of the Cu dyz 
orbital; likewise, correlation effects must involve an on- 
site repulsion Up rather than Urf and therefore should 
be smaller than might have been anticipated. Likewise, 
the exchange splitting of this band will be significantly 
different from that of the other d orbitals. 

The spin wave dispersion curves have been measured 
by Boehm et al. Their fit requires important nnn ex- 
change couplings, consistent with our finding that nnn 
hopping is essential to account for the dispersion. Al- 
though our hopping parameters are small, the exchange 
constants they obtain are very small, |J| <0.4 meV. Ap- 
plication of the superexchange expression J=4i^/U leads 
to unphysically large estimates of U«70-100 eV, indicat- 
ing that exchange coupling requires more careful consid- 
eration. 

Relative to the paramagnetic system, the FM state 
gains 70 meV/f.u. AF chains lead to a slight lowering 
by another 1.5 meV/f.u., but the phase with oppositely 
aligned FM chains is lower still (consistent with obser- 
vation), 80 meV/f.u. below the unpolarized state. The 
density of states for these phases are shown in Fig||, re- 
vealing an insulating gap of 0.32 eV. As usual, correlation 
effects will widen the gap, but we are not aware of any 




FIG. 4. Illustration of the overlap of two densities ICyzP 
centered on second neighbors Cu ions. Cu and O sites are 
denoted by empty and filled circles respectively. Model Gaus- 
sian dyz and pa orbitals were used for this plot. 

published data. 

The moment of 0.92 per Cu02 unit is roughly 60% 
on Cu and 40% on O ions. This transfer of magnetic mo- 
ment from a transition metal ion to a ligand ion (almost 
0.2 on each O) is to our knowledge the strongest yet 
reported in a transition metal oxide compound [Q. 

Previous reports of moments lie in the 0.02-0.10 /is 
range. Q To illustrate the importance of the O sites, we 
display the exchange potential V-\ — Vi for the ferromag- 
netic and antiferromagnetic chains in Fig. ^ There is a 
clear similarity, especially for the AF chain, to the T)*^ 
density in Figjj. Two other features should be noted: 
the exchange potential on O is comparable to that on 
Cu, and the exchange potential is of predominantly one 
sign for the entire V*^ orbital for both FM and AF or- 
dered chains. Unlike the strong dependence of O moment 
on the magnetic order, the size of the moment on the Cu 
ion itself (0.50-0.55 /is) is essentially independent of the 
overall magnetic ordering. 

Antialignmcnt of FM chains leads to the insulating 
band structure shown in Fig. 2(c). As expected, the 
two bands are described well by the eigenvalues of the 
system 

/il,l(fc) + iA ti,2(fc) \ 

where ti i = t2,2 contains all intrachain hopping, ti 2 = 
^2.1 contains the interchain hopping given above, and 
A=0.8 eV is the exchange splitting evident in Fig. 2(c). 
Since A is much greater than all of the hopping ampli- 
tudes and in fact similar to the PM bandwidth, the ex- 
change coupling between the antialigned chains results in 
very narrow band splits by A. 

Although the chains are nominally fully polarized, the 
coupling reduces the net moment to 0.92/iB/f.u. This to- 
tal moment is in excellent agreement with neutron scat- 
tering results, although the moment there was 
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FIG. 5. Exchange potential for ferromagnetic and antifer- 
romagnetic Cu-0 chains. Note for the AF case the close re- 
semblance to the T^y^ orbital density. The solid and dashed 
lines represent negative and positive values respectively. 

attributed solely to Cu. This is a remarkably large value 
for a quasi- ID chain, where quantum fluctuations 

should be large. The strong nnn coupling that we have 
uncovered account for the observation qualitatively: cou- 
pling between the V*^ effective orbitals, and therefore the 
spins, is really three dimensional, hence quantum fluctu- 
ation effect are vastly reduced. 

Recently two similar compounds with one-dimensional 
chains have been reported. Both Sro.73Cu02 and 
Cao.85Cu022 have Cu02 ribbons and order magnetically. 
[ p5[ Comparison with Li2Cu02 is difficult, however, be- 
cause the strongly differing doping level and the struc- 
tural disorder can lead to large changes in the magnetic 
behavior. Recent work on the Ca2+2:Y2_2;Cu50io 0] 
is suggestive, since susceptibility measurements are in- 
terpreted as indicating that each extra doped hole cre- 
ates a (non-magnetic) Zhang-Rice singlet. Our results for 
Li2Cu02 suggest that for the edge-sharing chains each 
extra hole occupies a V*^ orbital and therefore is weakly 
coupled to neighboring Cu ions. In this case the Cu ion 
and four neighboring po- orbitals would be non-magnetic, 
rather than having a Cu spin be compensated by neigh- 



boring O spins. Further experiments will be necessary to 
test our picture. 

To summarize, we have found that LDA provides a con- 
sistent picture of the magnetic properties and insulating 
character of the quasi-one-dimensional antiferromagnet 
Li2Cu02. Due to the formation of a strongly hybridized, 
and energetically isolated, combination of dyz and or- 
bitals, a large moment is transferred to the O ions. A 
simple single-band system results, but one in which sec- 
ond neighbors coupling exceeds nearest neighbors cou- 
pling and the electronic and magnetic behavior is three 
dimensional. 
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